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ABSTRACT 

We predict the amplitude of the gravitational redshift of galaxies in galaxy clusters using an N-body 
simulation of ACDM universe. We examine if it might be possible to detect the gravitational effect on the 
total redshift observed for galaxies. For clusters of mass M ~ 10 15 h _1 M©, the difference in gravitational 
redshift between the brightest galaxy and the rest of the cluster members is ~ 10 km s . The most 
efficient way to detect gravitational redshifts using information from galaxies only involves using the full 
gravitational redshift profile of clusters. Massive clusters, while having the largest gravitational redshift 
suffer from large galaxy peculiar velocities and substructure, which act as a source of noise. This and 
their low number density make it more reasonable to try averaging over many clusters and groups of 
relatively low mass. We examine publicly available data for 107 rich clusters from the ESO Nearby 
Abell Clusters Survey (ENACS), finding no evidence for gravitational redshifts. Test on our simulated 
clusters show that we need at least ~ 2500 clusters/groups with M > 5 x 10 13 h~ 1 M Q for a detection of 
gravitational redshifts at the 2a level. 

Subject headings: Cosmology: observations - large-scale structure of Universe 



1. INTRODUCTION 

General Relativity predicts redshift of photons due to a 
gravitational field. When a photon with wavelength A is 
emitted in a gravitational potential it will lose energy 
when it climbs up in the gravitational field and will con- 
sequently be rcdshiftcd. The redshift observed at infinity 
is given in the weak field limit by: 



AA 



(1) 



where AA, A$ are respectively the difference in wave- 
length, and difference in potential between where the pho- 
ton is emitted and where it is observed. 

If we consider galaxies as sources of the photons, the 
gravitational redshift effect is so tiny that we take it for 
granted that a measurement of the total galaxy redshift 
can be assumed to be the sum of Hubble expansion and 
peculiar velocities. In this paper we examine whether this 
is always the case, and in particular whether galaxies in 
galaxy clusters could have measurable values of z g . 

Since the gravitational potential depends on the mass 
distribution around galaxies, the gravitational redshift, if 
observable, should be most evident in dense environments. 
In an early study by Nottale (1976), the redshift difference 
between pairs of clusters was compared to the richness dif- 
ference. A supposed strong effect was found, with the pair 
member of higher richness having a systematically high 
redshift. However, when Rood & Struble (1982) rexam- 
ined this with a larger sample, their result showed no such 
correlation. Nottale (1990) discussed that the effect should 
be looked for in galaxies at the centers of galaxy clusters, 
by comparing their redshifts with those of galaxies at the 
cluster edges. Stiavelli & Setti (1993) carried out a re- 
lated test in individual elliptical galaxies, finding at 99.9% 
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confidence that elliptical galaxy cores are redshifted with 
respect to the galaxy outer regions, explaining this as a 
result of gravitational redshift. 

The study of gravitational redshifts in galaxy clusters 
was taken further by Cappi (1995), who modelled clusters 
using different density profiles including a de Vaucouleurs 
law. It was predicted that the gravitational redshift is non- 
negligible in very rich clusters. For example, the centers of 
clusters of masses 10 16 h _1 M Q should be redshifted by as 
much as 300 km s _1 with respect to infinity. Broadhurst 
& Scannapieco (2000) modelled the effect using a Navarro 
Frenk & White (1997) (hereafter NFW) density profile, 
and suggested that the gravitational redshift of metal lines 
in the cluster gas could eventually be used to map out the 
potential directly. As the gravitational redshift is sensitive 
to the distribution of mass in the innermost regions of clus- 
ters, it could be used as a probe to constrain the amount of 
dark matter there. Gravitational redshifts would provide 
complimentary information to gravitational lensing (e.g., 
Sand et al. 2003), as unlike lensing they do not depend on 
the mass density projected along the line of sight. 

Here we use an Y-body simulation made publically 
available by the Virgo Consortium (Frenk et al. 2000) 
to estimate the magnitude of the effect of gravitational 
redshifts on galaxy clusters in a ACDM universe. We ex- 
amine possible observational strategies and determine if 
galaxy gravitational redshifts could be detected with a rea- 
sonable number of clusters. By using the numerical simu- 
lation, we will be able to study the effect of substructure 
in the density and the potentially complex velocity field 
of realistic clusters. We will see if Nottale's suggestion of 
measuring the difference in gravitational redshift between 
the central galaxy and galaxies at the edge of the cluster is 
realizable in practice. The potential wells should be deeper 
for the most massive clusters, which means larger gravi- 
tational redshifts. However, massive clusters are rare, so 
that there may not be enough large mass clusters in reality, 
something which we will investigate. 

The gravitational redshift of the galaxy closest to the 
center of the potential is always positive with respect to 
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Fig. 1. — Some cluster examples. The top left panel shows the 
most massive cluster in the simulation (M = 1.4 X 10 15 h — 1 Mq); top 
right: M = 8.3 X lO^h^M©; bottom left: M = 1.1 X lO^h" 1 ]^; 
bottom right: Af = 1.2 X 10 15 h _1 M Q . Each panel is 8 /i _1 Mpc 
wide. 



the other galaxies, whereas the Hubble velocity component 
and peculiar velocity components can have either sign. We 
extract the gravitational rcdshift information by summing 
the total redshift from galaxies in many clusters so that 
the noise from Hubble expansion and peculiar velocities 
is averaged out. Small mass clusters may be better can- 
didates than massive ones in this sense because they are 
more abundant in the Universe. The effective mass range 
for detection is also discussed in this paper. 

The structure of the paper is as follows. In §2, the nu- 
merical simulation used as our model universe is described. 
In §3, we discuss the magnitude of the gravitational red- 
shift and the number clusters needed in order to detect 
the effect observationally. In §4, we briefly examine con- 
straints on gravitational rcdshifts from an observed sample 
of cluster galaxies, take from the ENACS survey. Finally, 
in §5, we summarize our results and conclude. 

2. SIMULATIONS 

2.1. Model Universe 

Our model universe is an output of an iV-body simula- 
tion run by the Virgo Consortium(Frenk et al. 2000), in 
which discrete particles represent dark matter. The back- 
ground cosmology is dominated by a cosmological constant 
(l^A = 0.7) and the matter content is 85% Cold Dark 
Matter (CDM). No gas dynamics is included and dissi- 
pationless particles are used to model the baryons also. 
The output we use is of the ACDM model at z = 
(O^m = 0.3, h = 0.7, as = 0.9). The box is periodic, of 
side length is 239.5 h~ 1 Mpc and contains 256 3 particles of 
mass 6.86 xlO 10 h^h^MQ. 

2.2. Cluster and Galaxy Selection 

For group finding, we use the friends-of-friends (fof) 
method (e.g., Huchra & Geller 1982). A particle is defined 



Fig. 2. — Slices through the gravitational potential profile of differ- 
ent mass clusters (same as Figure 1). Contours are of gravitational 
redshift with respect to oo and are labelled in km/s. 

to belong to a group if it is within some linking length (6) 
of any other particle in the group. We select clusters by 
using b = 0.2 where b is the linking length as a fraction 
of the mean particle separation. This definition of dark 
matter haloes was shown by Jenkins et al. (2001) to yield 
a mass function with a univeral form. Running the fof al- 
gorithm again, this time with b = 0.05, we find the largest 
group in each cluster, and define this to be the most mas- 
sive or central galaxy. When we refer to "galaxies" apart 
from the most massive galaxy later in the paper, we refer 
to particles randomly sampled from those cluster member 
particles which do not fall in the central galaxy halo. 

In Figure 1 we show some examples of rare rich clus- 
ters selected from the simulation in this way. The most 
massive cluster in the entire simulation has mass ~ 1.4 x 
10 15 h~ 1 MQ. Because the simulation includes non-linear 
clustering and clusters form from the bottom up, we expect 
each cluster to have its own substructure. Even though we 
are at z = and one might have expected some relaxation 
to have occurred, the substructure is very obvious. The 
cluster in the bottom right plainly consists of two clusters 
with similar masses merging. If such a merger was seen 
face on, then observationally it would be possible to sep- 
arate the clusters. However, if the clusters lie one behind 
the other along the line of sight, this would not be possible. 
In analyzing our simulated clusters, we are conservative 
and do not assume that clusters with much substructure 
could be removed, so that we use all clusters selected with 
our fof algorithm. 

2.3. Calculation of Gravitational Redshifts 

As a measure of gravitational redshift, we adopt velocity 
units, so that the quantity is given by: 

A$ 

Av g = Az g c= , (2) 
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where $ is the gravitational potential. For each particle 
in the simulation, we calculate the gravitational potential 
with respect to infinity using a tree algorithm (Barnes & 
Hut 1986). When doing this, we soften the gravitational 
potential in order to avoid noise from particles very close 
to one another, using 5 h _1 Kpc equivalent Plummer soft- 
ening (e.g., Aarseth 1963). 

In Figure 2 we plot contours of gravitational redshift in 
some examples of galaxy clusters (the same clusters whose 
density distributions were shown in Figure 1.) Each panel 
shows the two dimensional slice of gravitational potential 
in the xy plane that passes through the cluster center of 
mass. The substructure in the cluster density field is also 
somewhat reflected in the potential profile. However, the 
potentials are smoother, and even the merging clusters 
in the bottom right panel do not show strong evidence 
of two minima in the potential. The gravitational red- 
shifts, which are related to the potential will therefore vary 
smoothly across the cluster, without much substructure. 
The substructure in the density distribution will however 
cause noise in the Hubble flow part of the total redshift 
(see §3.2). 

In Figure 3 we show the effect of gravitational redshifts 
on the positions of particles when they are plotted in red- 
shift space. We choose y as the line of sight. In this figure 
we have not included the effect of peculiar velocities so 
that the much smaller gravitational redshift effect can be 
clearly seen. The value of Av g for the central galaxy is 
<~ AOkm/s, which causes it to be displaced 0.4 /i _1 Mpc 
further away from the observer than in real space. The 
gravitational redshift gets smaller rapidly as we move away 
from the center of the cluster and particles at the edges of 
the cluster mass distribution have Av g values wrt infinity 
of ~ 5 km s _1 . 

3. RESULTS 
3.1. Comparison with an Analytic Profile 
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Fig. 3. — The effect of gravitational redshift for the most massive 
cluster. The line of sight lies parallel to the y-axis. The lines start at 
the position of a galaxy before adding its gravitational redshift and 
end at the position including the gravitational redshift. Peculiar 
velocities are not added in the plot. 



We compare the gravitational redshift from the simula- 
tion with an analytic prediction using a spherically sym- 
metric density profile. We are interested in the difference 
in gravitational redshift between the central galaxy and 
the other cluster members and we will compare our ana- 
lytic results with those from the simulation. The density 
profile we adopt is given by (Hcrnquist 1990): 



p(r) oc 



1 



r(a + r) 3 



(3) 



where a ~ i? e /1.8153. The spatial gravitational redshift 
profile implied by this is (Cappi 1995): 



V H = 



GM 

c(r + a) 



(4) 



The gravitational redshift of the central galaxy or of the 
cluster is then given by integrating the mass weighted Vh 
value over the radius of the galaxy or the radius of the 
cluster, as in the following expression: 



P (r)V H {r)dV _ GMlo 9 



(a+r) 



■dr 



(a+r)' 



■dr ' 



(5) 



Io g P(r)dV c Q 

where here r g is the galaxy radius. 

We use the trapezoidal rule to integrate Eq. 5 numer- 
ically. For the characteristic radius of the central galax- 
ies and the clusters, we take the mass weighted mean ra- 
dius of galaxy (or cluster) member particles. With this 
definition, the average characteristic radius of clusters in 
the simulation with mass greater than 10 13 h~ 1 Mpc is 
0.338 ft- _1 Mpc, and the average radius of galaxies in 
those clusters is 0.084 /i _1 Mpc With a threshold in mass 
of 10 14 /i _1 Mpc the average radius of clusters and galax- 
ies respectively is 0.713 /i _1 Mpc and 0.154 /i _1 Mpc. The 
characteristic radius of the largest cluster in simulation is 
1.9 h^Mpc. 

For the gravitational redshift of the most massive galaxy, 
Eq. 5 has been integrated up to the radius of the most mas- 
sive galaxy, and for the gravitational redshift of the rest 
of members, the size of the most massive galaxy and the 
radius of the cluster are used as the lower and upper lim- 
its of integration respectively. The difference between the 
two values are calculated and the results from the simula- 
tion and from Eq. 5 are plotted in Figure 4. We can see 
from the figure that the values of central galaxy gravita- 
tional redshift from the simulation are broadly consistent 
with the estimate from analytic Hernquist density pro- 
file, although the estimate from analytic profile is slightly 
larger (by about 25%) than the result from the numer- 
ical simulation. We investigated if this could be due to 
the substructure of clusters, i.e., we checked whether Av g 
is large when Ay c . m . is small, which might be expected 
if concentricity of the galaxy and cluster yields a greater 
gravitational redshift. We find no such effect, however. 
A possible explanation of the difference could follow from 
the known fact that the density profiles of halos are better 
fitted by the NFW profile which gives p oc 1/r 3 for large r 
rather than the steeper (p oc 1/r 4 ) Hcrnquist profile. The 
latter would lead to a greater difference in the potential 
of the outer and central parts. In any case, the scatter we 
see in Figure 4 is quite small, indicating that substructure 
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Fig. 4. — Difference in gravitational redshift between the most 
massive galaxy and the rest of cluster members, x axis: from Hern- 
quist analytic profile using Eq 5, y axis: from simulation. The 
straight line is where Vg hernquist = Vg simu i ation . 



and deviations from spherical symmetry do not disrupt the 
agreement with the analytic results. Also apparent from 
the plot is the range of gravitational redshifts which will 
be relevant in our study from ~ 1 km s _1 to ~ 10 km s _1 . 

3.2. Gravitational Redshift in Redshift Space 

The total redshift velocity difference we observe between 
a galaxy and cluster is given by: 



Av t . 



Av g + AyxH + Av, 



(6) 



where Ay x Ho is the Hubble expansion part and Av is the 
peculiar velocity. For each component, we have calculated 
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Fig. 5. — The difference in gravitational redshift velocity (top 
left), center of mass (top right), and center of mass velocity (bot- 
tom left) of the central galaxy from the rest of the cluster members 
as a function of velocity dispersion. The bottom right panel shows 
the total redshift velocity, which is the sum of the 3 previous com- 
ponents. 1,000 randomly chosen clusters are shown. 



the difference in redshift between the most massive galaxy 
and the other cluster members (see also §3.3). There will 
also be measurement errors on the redshift, which we shall 
cover later. In Figure 5 we plot each redshift component 
separately and the total redshift as a function of the cluster 
line-of-sight velocity dispersion. In order that these values 
may be related to cluster mass, we have fitted a power 
law to the relationship between cluster mass and velocity 
dispersion, finding it to be: 



logM = (2.56 ±0.04) logo- + 7.12 ±0.08, 



(7) 



where M is in units of h _1 MQ and the units of a are 
km s _1 . The slope of the relation is not too far from the 
value of 2 expected for virialized objects. The massive, 
high velocity dispersion clusters have the largest central 
galaxy gravitational redshifts. The rest of the redshift 
components exhibit larger scatter as velocity dispersion 
increases. 

The cluster in the bottom right panel in Figure 1 has 
more substructure than the others and this is reflected 
in the large redshift difference between the most massive 
galaxy and the center of mass of rest of members. The 
position differences for all the panels are Ay c , m , — 0.0698, 
0.049, -0.181 and 1.309 /i _1 Mpc respectively (clockwise 
from the top left in Figure 1). Compared to the veloc- 
ity difference, however, this has a small effect on the total 
redshift (see top right panel in Figure 5). Figure 5 shows 
that most scatter comes from the velocity difference be- 
tween the central galaxy and the rest of the cluster. 

The plot shows how small the gravitational redshift is 
compared to the other redshift components. In the bottom 
right panel of Figure 5 where we show the total redshift, 
it is obvious that we cannot see the effect of gravitational 
redshifts by eye on individual clusters or overall in a scatter 
plot such as this one. 

The redshift due to the position and velocity differences 
will be the most important source of noise when we try 
to extract the gravitational redshift. If the most massive 
galaxy were lying at rest at the centroid of the cluster, then 
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Fig. 6. — Difference in redshift between central galaxy and cluster 
centroid as a function of cluster velocity dispersion a. We show 
the scatter in center of mass velocity, Hubble component, and total 
redshift, as in Figure 5. The lines show the values which enclose 68% 
of the distribution and symbols the Icr standard deviations. Error 
bars (Poisson errors) are shown for the center of mass velocity. 
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the accuracy of our measurement of the gravitational red- 
shift difference between it and the other cluster members 
would be limited only by Poisson noise in computation 
of the mean redshift of the cluster. Because the central 
galaxy is usually moving with respect to the mean of the 
other cluster members, and is usually not at the cluster 
center of mass, we have this additional source of noise, 
shown in the bottom left panel of Figure 5. Averaging 
results over many clusters in order to reduce the error on 
the mean v g is the approach we will take in this paper. 
An interesting question is how best to do this. For the 
clusters with the greatest mass, v g will be large, but so 
will the noise per cluster. On the other hand, v g for small 
mass clusters is small but the peculiar velocity and posi- 
tion noise is small as well. We will discuss this more later 
in the paper and determine the optimal range of cluster 
mass for detection of gravitational redshifts. 

One can ask how one might expect the dominant noise 
component AV c . m . to scale with the cluster velocity dis- 
persion a. If the central galaxy were a random galaxy in 
the cluster, then the scatter in AV c . m . should be approx- 
imately a. However, we expect it to be substantially less 
than this, as the central galaxy, while not at rest is still 
moving less with respect to the cluster center of mass than 
a random galaxy. In order to quantify this, we have plot- 
ted in Figure 6 the standard deviation of AV c . m . against a. 
We find standard deviation of AV c . m . to be approximately 
0.35(T. We have also plotted the other source of noise from 
the Hubble component Ay c . m . x H , as well as the total, 
AVtot- Interestingly, the total noise on the central galaxy 
position, AVtot is slightly lower than AV c . m .. On inves- 
tigation, we have found that this is because of residual 
coherent infall of particles into the cluster, so that clusters 
on the far side, with +ve Ay cm . x Hq have -ve values of 
AV c . m . an d vice versa. In Figure 6, we also compare the 
directly computed 68% interval to the standard deviation, 
finding that the distribution of AV to t is close to Gaussian. 

3.3. Estimate of the Number of Clusters Needed to 
Observe Gravitational Redshifts 

One simple method for detection of gravitational red- 
shifts would involve finding the redshift difference between 
the central galaxy and the other cluster members, after av- 
eraging over a large number of clusters. We have seen in 
the previous section that the noise from substructure is 
small, and that this should be feasible given enough clus- 
ters. In this subsection we investigate how well this could 
be done in practice. In the next subsection, we will include 
information on all galaxies in the cluster rather than just 
splitting the galaxy population into "central galaxy" and 
"other". 

With the difference in gravitational redshift computed 
from the simulation, we make an estimate of the number 
of clusters needed for a detection. We try two cases, ran- 
domly selecting either N ga i — 50 galaxies per cluster or 
N g ai = one tenth of the number of dark matter particles, 
so that N ga i — M c i uster /7 x 10 11 h _1 M Q . We assume that 
the noise from the difference in position and velocities con- 
tributes independently. Including the measurement uncer- 
tainty A meas ., we have 
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And the error on the mean is given by 



(9) 
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We have split the clusters in the simulation into bins 
by mass, and N c here is the number of clusters in the 
simulation in a particular bin. For each bin we com- 
pute the number of clusters (N ) actually needed for a 
detection of the gravitational redshift, with a given sta- 
tistical significance. For example, for clusters in the mass 



Otot 



bin centered on ~ 8 x 10 14 h 1 M P) , a <~ 300 km s 



and 



Av a 



7 km s 



10 4 clusters we expect our cr- 
which is a ~ 2a detection of 



With ~ 

ror bar to be <~ 3 km s _1 
gravitational redshift. 

We might expect that the number of clusters should 
decrease for high mass because the gravitational redshift 
is proportional to the cluster velocity dispersion. How- 
ever, for these clusters, the error bars also increase be- 
cause the dispersion in velocity difference between central 
galaxy and others is getting larger. Figure 7 shows how 
many clusters we need to detect the gravitational redshift 
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Fig. 7. — The number of clusters needed to detect graviational 
redshift as a function of cluster mass. We show results using either 
50 randomly selected galaxies per cluster (top) or using a number 
of galaxies proportional to the cluster mass (bottom), both for a 2a 
level detection. For the latter case, the number of galaxies used is 
equal to 0.1 times the number of dark matter particles in it, so that 
for a IO^Ii^Mq we use 15 galaxies. Measurement uncertainties 
of 30 km/s, 100 km/s and 300 km/s are included. The thick line 
represents a 4cr detection with no measurement uncertainty. Note 
that the graph flattens for large cluster masses. 
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at the 2<T and 4cr levels as a function of cluster mass. We 
choose the measurement uncertainty to be 30 km/s, 100 
km/s and 300 km/s as best, typical and worst cases (e.g., 
Stoughton et al. 2002). Given relatively small measure- 
ment uncertainty, the number of clusters needed decreases 
at first and then stays roughly the same (~ 10 4 clusters for 
A meas . = 30 km/s) after M > lO^h" 1 ]^©. Since we use 
average values, it does not make much difference whether 
we use 50 galaxies per cluster or a number proportional to 
the cluster mass. 

In order to understand the behaviour seen Fig 7, we can 
ask what should happen if the redshift noise on the cen- 
tral galaxy from peculiar velocities and position differences 
were proportional to the cluster velocity dispersion (as- 
suming we are averaging over a fixed number of galaxies in 
each cluster). As the gravitational redshift is proportional 
to cluster velocity dispersion also (Figure 5) , then the num- 
ber of clusters for a detection of a given significance should 
be constant. This is approximately what is seen for clus- 
ters of mass M greater than ~ 10 14 h~ 1 M Q . For smaller 
clusters, the noise from the displacement of the central 
galaxy (the Hubble velocity component, Ay c . TO . x Hq) is 
becoming a non negligible fraction of the total noise, and 
its effect means that more clusters must be averaged over. 

Note that, as mentioned above, only the mass-averaged 
gravitational redshift has been used here and it has not 
been taken into account that the gravitational redshift has 
a gradual radial dependence. More information is there- 
fore available which could in principle aid a detection. We 
discuss this in the following subsection. 

3.4. Gravitational Redshift Profile 

We now calculate the profile of gravitational redshifts 
averaged over many clusters. As in the previous section, 
we take the central galaxy to be our zero point. The other 
galaxies in the cluster are binned as a function of their im- 
pact parameter from this galaxy (r = \J Ax 2 + Az 2 , where 
y is the line of sight) In each bin, and for all clusters, we 
have summed the gravitational redshift difference between 
the galaxy and the central cluster galaxy. We then divide 
by the number of galaxies in each bin (as mentioned pre- 
viously the "galaxies" which are not the central one are 
actually particles). 

Figure 8 shows the resulting gravitational redshift pro- 
file for clusters within different mass ranges as a function of 
the impact parameter r. All of the clusters in the simula- 
tion have been used, and it should be noted that each mass 
bin does not an have equal number of clusters contribut- 
ing to it. Also, at large radii, within each mass range, only 
the more massive clusters contribute, because the smaller 
clusters do not have any galaxies out that far. As ex- 
pected, gravitational redshifts increase for large mass clus- 
ters. This result is consistent with that of Cappi (1995) 
who used analytic profiles, in that the function decreases 
smoothly with a maximum at the center. 

As we have described in section 3.1 our central galaxies 
have a mean radius of ~ 0.1 /i _1 Mpc, which means that 
within the central few bins, we do not expect the gravi- 
tational redshifts of the galaxy to rise very steeply. If the 
central galaxy were modelled as a point (as was done by 
Broadhurst & Scannapieco 2000), the profile could become 
steeper. The dark matter halos in this simulation will have 
profiles consistent with the NFW form, so that we do not 



expect the central profile to be very steep even if the cen- 
tral galaxy is allowed to be much smaller, at least in this 
dark matter only run. The gravitational redshift profile 
for NFW clusters is shown in Broadhurst & Scannapieco 
2000., and flattens off towards the center. According to El- 
Zant ct al. (2003), any baryonic dissipational effects can 
cause either a steeping or flattening of the density pro- 
file towards the center, which will increase or decrease the 
gravitational redshift, but this is ignored in the present 
study. 

In order to use the redshift profile information to make 
a detection of gravitational redshift, we must average over 
a large number of clusters to reduce the noise from pecu- 
liar velocities and Hubble expansion. As before, we add 
the total redshifts for many clusters so that the redshifts 
from position and velocity are averaged out, leaving only 
the gravitational redshift. This time, this will be done 
in bins of impact parameter. We can then estimate how 
many clusters we need to observe in order to ascertain to 
a given significance level that the gravitational redshift is 
inconsistent with zero. 

First, given a sample of simulated clusters in a certain 
mass range, we calculate the gravitational redshift profile 
in bins of impact parameter and the matrix of covariance 
between bins. The technique of Singular Value Decompo- 
sition (Press et al. 1992) is used to invert the matrix, and 
then the value of \ 2 is computed: 

^^(xi-x^ixj-Xj) (12) 

AT 

where Xi is the gravitational redshift, dj is covariance ma- 
trix. As we are interested in the possibility of detection, 
we set Xi = Xj = exactly so that \ 2 is t nen the differ- 
ence A% 2 from a model with zero gravitational redshift. 
Because dj is proportional to 1/ ^NiNj, where Ni is the 
number of clusters contributing to r bin i, \ 2 is propor- 
tional to the total number of clusters, if Ni — Nj. The 
number of r bins which could be used depends on the an- 
gular extent of the cluster. Our criterion is to truncate at 
the outermost bin with less than 20 clusters contributing 
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Fig. 8. — The gravitational redshift profile averaged over many 
clusters as a function of impact parameter from the most massive 
galaxy. The gravitational redshift of the central galaxy is subtracted 
from the average gravitational redshift of the rest of cluster mem- 
bers. 



7 



10" 



10 



° 10' 



10' 
10 




10" 
mass (h _1 M G ) 



10" 



10' 



10 



10' 
10 




10" 
mass (h _1 M G ) 



10" 



Fig. 9. — top: Number of clusters needed to detect gravitational 
rcdshifts at the 2<r (circle) and 4<r ( X ) levels as a function of cluster 
mass. The x 2 has been calculated with (dotted line) and with- 
out the (solid line) covariance matrix. We use all particles in each 
cluster, bottom: Number of clusters needed to detect gravitational 
redshifts as a function of mass when we are restricted to 100 galax- 
ies per cluster. Note that the graph flattens for the mass above 
~ 3 X 10 1A h^ 1 Mq when only diagonal terms are used. Measure- 
ment uncertainties are not included. 



and not use the larger r bins to build the covariance ma- 
trix. Additionally, when the number of clusters in a given 
mass range is small, the matrix is too noisy to invert. We 
consider the results invalid when x 2 calculated with the 
covariance matrix is greater than x 2 computed using only 
the diagonal elements of the covariance matrix. This oc- 
curs for the clusters with M > 5 x 10 14 h~ 1 M Q since there 
are only 4 in the simulation, and for this mass range we 
use the diagonal elements of the matrix only. 

The top panel in Figure 9 shows the number of clusters 
needed for detections of non-zero gravitational redshifts 
at the 2ct and 4<r levels as a function of mass. No mea- 
surement uncertainties are included in this plot, x 2 values 
with only the diagonal elements were also calculated and 
those results are shown as additional lines in the plot. In- 
cluding the proper covariances between bins reduces the 
sensitivity of the test, so that approximately a factor 2 
more clusters are needed. The plot has slightly different 
behaviour to Figure 7, as the number of clusters needed 
continues to decrease for high masses. This is basically be- 
cause we have used all galaxy information in each cluster. 
This is obviously not very realistic, as up to 20000 galax- 
ies (particles) have been averaged over in the most massive 
clusters. We have therefore recomputed our results assum- 
ing that only 100 galaxies per cluster are available. 

When we do this as shown in the bottom panel in Fig- 
ure 9, we find that the curve flattens for high mass clusters 



and even goes up slightly when the effect of covariances 
is included. This means that observing clusters of mass 
~ 5x 10 13 h _1 M Q seems to be the most efficient way to 
detect gravitational redshifts, at least when using galax- 
ies as tracers of the potential. For higher mass clusters, 
we have larger gravitational redshift effects but also have 
more noise. 

Note that when only the diagonal terms in the covari- 
ance matrix are used in % 2 calculation the flattening of 
the graph is consistent with the simpler estimator in §3.2 
(Figure 7). However, both the mass threshold and the 
number of clusters needed to detect gravitational redshifts 
are smaller here, by about a factor of 3. This is because 
we have more information in the sense that we use the 
gravitational redshift profile in many r bins while in §3.2, 
there are effectively only two bins - the central galaxy and 
the rest of the galaxies in each cluster. 

For a more direct comparison with Figure 7, we then 
use 50 galaxies per cluster and one-tenth the number of 
particles while including measurement uncertainties in the 
covariance matrix (Figure 10). When the uncertainty is 
small, the result is consistent with Figure 7 in that with 
fixed number of galaxies per cluster, the curve flattens 
for high mass clusters. However, when the uncertainty is 
large, it appears that the measurement uncertainty is the 
dominant factor deciding the number of clusters needed 
regardless of the number of galaxies per cluster. 

We have seen that for a 2a detection of gravitational 
rcdshifts we expect to need a sample of ~ 5000 galaxy 
clusters/groups with 50 measured redshifts per cluster ( 
with measurement error A meas . = 30 km/s). A statisti- 
cally correct detection would involve making use of the 
covariance matrix of the noise between bins, which could 
be done from our model clusters, or else computed directly 
from the observational sample itself. This required num- 
ber of clusters/groups is large, but may be just within the 
range of present day redshift surveys. We will return to 
this in our discussion (§5), as well as exploring results from 
a recent survey below. Alternatively, as stressed by Cappi 
(1995), the observed absence of an effect could be used to 
place limits on the gravitating mass of clusters. 

4. RESULTS FROM THE ENACS SURVEY 

The ESO Nearby Abcll Cluster Survey (ENACS, Kat- 
gert et al. 1998) is a redshift survey of rich Abell clusters 
in the southern sky. The catalog contains redshift, mag- 
nitude and position information for 5634 galaxies in 107 
clusters. The mean velocity dispersion of these clusters is 
<~ 700 km s" 1 (Biviano et al 1996), so that based on the 
work in §3 we only expect there to be a small (<~ 5 km s _1 ) 
redshift difference on average between the central galaxy 
and the other cluster members. We will compute this 
statistic in order to check consistency with our expecta- 
tions. 

We take the ENACS data and find cluster centers by 
placing a cylinder of length 3000 km s _1 along the line of 
sight and radius 2.0 /i _1 Mpc around each galaxy. We find 
the center of mass of each cylinder (weighting each galaxy 
by the inverse of the selection function), recenter it on this 
point and iterate until convergence. We eliminate all over- 
lapping cylinders, keeping those with most galaxies. We 
then rank all cylinders by the weighted number of galax- 
ies and end up with a list of 107 rich clusters and their 
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Fig. 10. — top: Number of clusters needed to detect gravita- 
tional rcdshifts at the 2a level as a function of cluster mass. The 
measurement uncertainties of 30 km/s, 100 km/s and 300 km/s are 
included. We use 50 particles per cluster (top) and tenth the size of 
cluster (bottom). Note that the graph flattens for the mass above 
~ 3 X 10 13 /i~ 1 Mq. The thick line represents the 4<r level with no 
measurement uncertainty. 



members. 

Because we do not expect to be able to detect a gravita- 
tional redshift signal, we only carry out the simplest test, 
measuring the difference in redshift between the bright- 
er ^ 
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Fig. 11. — Difference in redshift between the brightest cluster 
galaxy and the other cluster members, for 107 rich Abell clusters 
from the ENACS. A Gaussian curve with the same standard devia- 
tion, a = 486 km s _1 is also shown. 



est galaxy in each cluster and the other cluster members. 
The results are shown as a histogram in Figure 11. The 
rms difference between the brightest and other cluster 
galaxies is 486 ± 47 km s _1 (Poisson error). In the sim- 
ulations, for a sample of clusters with the same mean ve- 
locity dispersion we expect an rms difference of 255 km s _1 
(from Hubble flow and peculiar velocities.) This is rather 
lower than the observed value, and in fact it takes a simu- 
lated sample with a mean velocity dispersion of 970 km s _1 
to have the same rms difference (by applying a lower 
mass threshold of 9.3 x lO^h^M©. This could be a 
sign that there may some differences between the struc- 
ture of observed and simulated clusters. It is likely that 
this is at least partly due to observational measurement 
errors, as the quoted uncertainties of ~ 100 km s , when 
added in quadrature would increase the rms difference to 
275 km s -1 . The mean of the histogram for the ENACS 
clusters is —66 km s _1 (a blue shift), with a la uncertainty 
of 47 km s _1 so that as expected there is no detectable ev- 
idence for gravitational rcdshifts from this small sample. 

We note that Cappi (1995) also carried out a similar 
same test with a sample of 42 clusters with cD galaxies, 
and found a small blue shift, —39 ± 32 km/s. Although 
in the ENACS sample we have simply measured the red- 
shift difference of the brightest galaxy without regard to 
whether it is a true cD or not, we find essentially consistent 
results. On the other hand, Broadhurst and Scannapieco 
(2000) find a mean cD redshift of 260 ± 58 km s _1 for a 
sample of 8 clusters, of which 5 have cooling flows. This 
value is an order of magnitude larger than what we expect 
from our simulations. These authors have pointed out that 
the central potential could have been made steeper over 
time by cooling flows. 

5. SUMMARY AND DISCUSSION 

We have explored the gravitational rcdshifts of galaxies 
in the potential wells of simulated galaxy clusters. In par- 
ticular we have concentrated on the possibility of detect- 
ing gravitational redshifts from surveys of cluster galaxies. 
Our main findings are as follows: 

• For clusters of <~ 10 15 h _1 M Q , the difference in red- 
shift between the central galaxy and the other cluster 
members is ~ 10 km s _1 . 

• The gravitational potential of clusters does not ex- 
hibit much substructure, even when the clusters are 
not relaxed. The most important source of noise on 
a measurement is from galaxy peculiar velocities. 

• The most sensitive method for detecting gravita- 
tional rcdshifts with galaxy data is to use the grav- 
itational redshift profile as a function of impact pa- 
rameter. 

• For a given number of galaxy redshifts measured, 
the most efficient strategy would be to target clus- 
ters of <~ 5 x 10 13 1i -1 Mq, which have lower noise 
from galaxy peculiar velocities than more massive 
clusters. 

• For a 2a detection, we require ~ 5000 clusters with 
M > 5 x 10 13 h _1 MQ when the measurement un- 
certainty is 30 km/s and ~ 2500 clusters when the 
measurement uncertainty is negligible. 
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Our result for the shape of the gravitational rcdshift 
profile are consistent with that of Cappi(1995) and Broad- 
hurst and Scannapieco (200). However, the masses of clus- 
ters in our simulation are generally smaller than those used 
as examples by these authors. For example Cappi (1995) 
found that gravitational redshift with respect to infinity of 
the central regions of a cluster of mass M = 5x lO 15 h _1 M0 
to be > 100 km s _1 and predicted the effect for even more 
massive clusters with M = 10 16 h _1 M s (v g ~ 300 km s" 1 ). 
Unfortunately, even though the gravitational redshift in a 
cluster with mass 10 16 h _1 M Q is large, there is little chance 
of finding such a massive cluster in the Universe. 

We can roughly estimate the number of clusters in the 
universe whose mass is M = 8 x 10 14 h _1 M Q or greater. 
First, we approximate the mass function of clusters given 
by (Bahcall & Bode, 2003) as a power law: 

logn(> Mi. 5)Com = 8 x lO^h^M©) w -2z - 6.75, (13) 

The comoving distance as a function of redshift is given 
by (e.g., Hogg 1999): 

f z dz' 

D - D -L W) (14) 

where E(z) = \/fi M (l + z) 3 + fi fc (l + z) 2 + S1 A and D H = 
3000 h^Mpc. We take n M = 0.3 and fi A = 0.7. We then 
find the total number of clusters in the universe (up to 
z = 1, although raising the z limit makes essentially no 
difference) : 

(15) 

Integrating Eq. 15 numerically we find that there are only 
~ 600 clusters (M > 8 x lO^h^M©) observable in the 
entire universe. We have seen, however in §3, that in any 
case the larger noise from peculiar velocities in more mas- 
sive clusters makes a detection of gravitational redshifts 
problematic anyway. It might be better to use lower mass 
clusters in order to detect the gravitational redshift, even 
though effect is small. 

As for observed samples of clusters with lower masses, 
the situation is more promising. For example, the Sloan 
Digital Sky Survey has been used to find ~ 1,000 clus- 
ters in total with z < 0.15 and M > 5 x 10 13 h _1 M Q so 
far (Chris Miller, 2003, private communication). For a 
more on cluster finding in the SDSS, see Nichol (2003) 
and Annis et al. (2003). Published samples of SDSS clus- 
ters include those of Goto et al. (2002) and Bahcall et al. 
(2003b). Merchan and Zandivarcz (2002) have constructed 
a galaxy group catalog from the 100K data release of the 
the 2dF galaxy redshift survey (refs). They find <~ 1000 
clusters/groups with masses M > 5 x 10 13 h~ 1 M Q . The 
final 2dF data release contains ~ 2 x 10 5 galaxies (Col- 
less et al. 2003), which should enable construction of a 
group/cluster catalog with ~ 2500 members above this 
lower mass limit. As we have seen in §3.4, this may be 
close to the level required in order to make a 2a detection 
of the gravitational redshift. A catalog of similar groups 
selected from the final SDSS survey might allow a 4 a 
detection. 



Because the effect we are searching for is very small (a 
rcdshift on the order of a few km/s), we must be careful to 
consider possible systematic effects. The effects of pecu- 
liar velocities and cluster substructure have been included 
in our simulations, and averaging over large numbers of 
clusters produces reasonable results. One effect we have 
not included is full selection of clusters in rcdshift space. 
In the simulation, we simply define a cluster in real space 
with a linking length of our choice. In a survey, we select 
the clusters in redshift space and we include some contam- 
inating galaxies behind and in front of each cluster that 
are not members. This would act to increase the noise, 
but not by much, as they will be heavily outnumbered by 
cluster members. 

A real problem would be any effect which gives a system- 
atically positive or negative redshift for the central galaxy 
with respect to the other galaxies. If we select cluster 
galaxies using a cone of constant angular size rather than 
a cylinder of fixed radius, then more contaminating galax- 
ies will come from the larger volume behind the cluster 
than in front, and the average redshift of the galaxies apart 
from the central one will change. We have calculated the 
potential effect of this bias using the cluster-galaxy cross- 
correlation function data of Croft et al. (1999) to find 
the number of contaminating galaxies. We find a rela- 
tive blucshift of ~ 1 km s _1 for the central galaxy. This 
can be avoided simply by using a cylinder to select cluster 
members, however. 

An additional potential problem which goes in the oppo- 
site direction is the fact that with a magnitude limited sur- 
vey, the contaminating galaxies behind the cluster should 
be fewer in number because they are at a greater distance. 
The size of this effect depends on the steepness of the lumi- 
nosity function of galaxies. It could also however easily be 
eliminated by applying a local volume limit in the vicinity 
of each cluster, rejecting galaxies which would be too faint 
when placed at the far edge of the cluster. 

In this paper we have seen that realistically, gravita- 
tional redshifts are an extremely small effect, and that 
looking at the most massive galaxy clusters may not help 
in detection because the noise is much larger than for small 
clusters. If a detection is to be made using a galaxy survey, 
it will be important to carry out cross checks. For exam- 
ple, enough data needs to be a available that clusters can 
be divided by mass into more than one bin, to make sure 
that the effect is largest for the high mass clusters (even 
if the error bar is large). The more tracers of the clus- 
ters potential that are available, the better. Broadhurst & 
Scannapieco (2000) have shown that X-ray emission from 
intracluster gas could be used. Other probes of the poten- 
tial could include using redshifts of intracluster planetary 
nebulae (e.g., Fcldmeier et al. 2003). 

It may be marginally possible to detect the gravitational 
rcdshift of cluster galaxies with the data available today. 
However, more surveys of greater scope have been imag- 
ined for the future. Someday this will provide us enough 
information to separate the gravitational rcdshift effect 
from total redshift. 

We thank the anonymous referee for useful suggestions. 
RACC acknowledges support from the NASA-LTSA pro- 
gram, contract NAG5-11634. 
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